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CN ' ABSTRACT 

i Two approximations, namely the subseismic approximation and the anelastic approx- 

imation, are presently used to filter out the acoustic modes when computing low 
frequency modes of a star (gravity modes or inertial modes). In a precedent paper 
(Dintrans & Rieutord 2001), we observed that the anelastic approximation gave eigen- 
frequencies much closer to the exact ones than the subseismic approximation. Here, 
we try to clarify this behaviour and show that it is due to the different physical ap- 
. proach taken by each approximation: On the one hand, the subseismic approximation 

' considers the low frequency part of the spectrum of (say) gravity modes and turns out 

£S) \ to be valid only in the central region of a star; on the other hand, the anelastic ap- 

proximation considers the Brunt- Vaisala frequency as asymptotically small and makes 
no assumption on the order of the modes. Both approximations fail to describe the 
modes in the surface layers but eigenmodes issued from the anelastic approximation 
are closer to those including acoustic effects than their subseismic equivalent. 

We conclude that, as far as stellar eigenvalue problems are concerned, the anelas- 
+3 . tic approximation is better suited for simplifying the eigenvalue problem when low- 

frequency modes of a star are considered, while the subseismic approximation is a 
. . | useful concept when analytic solutions of high order low-frequency modes are needed 

J> . in the central region of a star. 
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1 INTRODUCTION 

When considering the low-frequency modes of a star, namely 
gravity modes or inertial modes, the compressibility of the 
fluid is often a side effect in the determination of eigenfre- 
quencies and eigenmodes; in other words, the dynamics of 
these modes may be simplified by neglecting the elasticity 
of the fluid or, equivalently, by filtering out acoustic modes. 
This is the aim of the subseismic and anelastic approxima- 
tions; the resulting equations for eigenmodes are much sim- 
pler than the original ones and very useful when dealing with 
the low frequency oscillations of rotating stars (e.g. Dintrans 
& Rieutord 2000). 

Recently, we compared these two approximations (Din- 
trans & Rieutord (2001) referred to as paper I hereafter). We 
found that in the two cases which we analysed, namely two 
polytropes, the anelastic approximation performed much 
better than the subseismic approximation. We attributed 
this behaviour to an inconsistency of the subseismic approx- 
imation but our argument turns out to be not general and 



Smeyers (2001) showed that, for low-frequency high order 
modes, the subseismic approximation gives the first order 
equations in regions not close to the surface of the star. 
These results prompted us to re-examine this question in or- 
der to clarify the origin of the different behaviour of the two 
approximations. For this purpose we will focus, in section 2, 
on two asymptotic developments: a first one where we use, 
as Smeyers (2001), the frequency as a small parameter and 
a second one where we use the Brunt- Vaisala frequency as 
the small parameter. These asymptotic developments will 
prove to be at the origin of each of these approximations 
and will allow us to clarify the physics attached to each of 
them. In section 3, using the same examples as in paper I, we 
will compare the approximate eigenfunctions to their exact 
counterparts and show the better behaviour of the anelastic 
approximation. Our conclusions are drawn in section 4. 
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2 THE ASYMPTOTIC EQUATIONS 

As was shown in I, both approximations imply Cowling's ap- 
proximation; we shall therefore neglect the perturbation of 
the gravitational potential and will start from the following 
equations: 



p' + div(pC) = 0, 



, p'\ N'^ 
w (, = V — - — SPe r , 



pg 



8P = c 2 Sp, 



(1) 



(2) 



(3) 



where we assumed a time-dependence of the form exp(itirf) 
and considered adiabatic oscillations. £ is the displacement; 
P' and p' respectively denote the Eulerian fluctuations of 
pressure and density whereas SP, 8p are their Lagrangian 
counterparts; thus we have 



S P = P +^r, 



with a pressure gradient satisfying the hydrostatic equi- 
librium dP/dr = — pg. Also, p is the equilibrium density, 
g = — ge r the gravity and 7 = (<9 In P/d In p)$ the first adi- 
abatic exponent. Finally, c 2 and JV 2 respectively denote the 
squares of sound speed velocity and Brunt- Vaisala frequency 
such as 



2 P 
c = 7 _, 



1 7 dr 



d\np 
dr 



(4) 
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However, if we use the expansion (5) into (7-8), we recover 
(6); therefore, (7-8) contain terms of higher order than (6). 

In fact, our expansion (5) breaks down near the origin 
r = where the regularity of the solutions (i.e. that £ oc 
r 1 ^ 1 ) is not insured. This comes from the fact that terms 
like f/r are no longer negligible compared to derivatives 
df/dr. 

This difficulty is avoided by Smeyers (2001) with the 
use of the variable r = ^^ l+v> f„ r ^ dr' instead of the 

s JO r' 

radial variable r. r serves as a fast variable while r, the slow 
variable, is assumed small compared to the scale of variation 
of the Brunt- Vaisala or the background density. Using this 
transformation, Smeyers (2001) shows that (8) is verified 
by the solution at leading order while (7) is approximately 
verified. 

Near the surface layers, Smeyers (2001) has shown that 
within this development, the Eulerian pressure perturbation 
is no longer negligible and that the subseismic approxima- 
tion does not apply. 

Thus the subseismic equations govern the oscillations 
of high radial order gravity modes in the central parts of 
the star. No constraint is imposed to the Brunt- Vaisala fre- 
quency and the equations to be solved are: 



div£=^£r, 

cr 



u£,= V ( — ) +N\ r e r 



2.2 The anelastic view 



2.1 The subseismic view 

As a first exercise we derive the equations verified by the 
low frequency gravity modes. 

We therefore assume that the frequency reads uj = eo>i 
and that ^ scales as e _1 with e <§C 1 since we focus on 
high radial order modes. Developing the dependent variables 
generically as 



(5) 



and using the classical expansion of the variables on the 
spherical harmonics, 



+00 +t 



£=0 m=-£ 

and dropping (£, m)-indices, we find that Pq = Pi = 0, 
Po = 0, Co = and that 



dr 



(rxo) 



N 



(6) 



^(r\i)=£(£+l)rxo. 



A system which is slightly different from the one ob- 
tained when setting P' = which yields the subseismic 
equations: 



Let us now turn to the anelastic approximation. In this case, 
it is more convenient to write the equations of motion for 
the velocity field v rather than the displacement £, that is, 



1 d 



iup' + 4-f (r 2 pu) -£(£ + 1)?!= 0, 



r 2 dr 



d P' 



dr \ p 
P' 



N 2 ( 
iujpg 



jP' - pgu), 



(9) 



tujrv = — - 



iui(P' — c 2 p) 



pc 2 N 2 



-u, 



where we used the following spherical harmonics decomposi- 
tion for v (for clarity, we still dropped in the previous system 
the £, m indices) 



+00 +1 

v(r,e,cj>) = J2 uUr)Y e m (0,Wr + vUr)VYr. 
1=0 m=-e 

We now assume that the Brunt- Vaisala frequency is 
vanishingly small; note that as this quantity often diverges 
at the surface of the star models, it is more appropriate to 
assume that ujn, the frequency of the lowest order gravity 
mode, is vanishingly small. Thus we write 

N(r) = en(r), w = sun, P' = Pq + eP[^ , ... 



dr 



Note that we make no assumption concerning the scale of 
(7) the perturbations which may be of order unity. First orders 
yield the equations 
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Figure 1. Normalized eigenfunctions £ with 1 = 2 and k = 5 for the homogeneous polytrope. The solid line shows the exact solution 
(w 2 ~ —3.61 X 10 -2 ) while the dashed and dotted lines correspond to its subseismic and anelastic approximations, respectively (with 
ojJ uba ~ —3.26 X 10~ 2 and i^ 2 n(!l ~ —3.53 X 10~ 2 ). On right, the surface layers have been magnified. 



Po = Po = 0, 
V • puo = 0, 

d(P{/p) 



ILOlUu = — - 
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«0, 



2 / \ 
C Pi 



2 2 

pn c 



-Uo- 



(10) 
(11) 
(12) 
(13) 
(14) 



from which we write the anelastic system 

div(p£) = 0, uj 2 i = V(P'/p) + N 2 £ r e r 

As for the subseismic approximation the perturbed equa- 
tion of state is eliminated; but from (14), we note that the 
Eulerian pressure perturbation is of the same order as c p' . 

The subseismic equations can be obtained by just drop- 
ping out the Eulerian fluctuation P' in SP, in the original 
equations (1-3); from (14), we see that this is not the case 
in the anelastic approximation. On the other hand, the fluc- 
tuation of density in the mass conservation equation can be 
neglected. 

We therefore see that the anelastic approximation ap- 
plies when the Brunt- Vaisala frequency is small compared to 
the acoustic frequencies but does not impose any constraint 
on the scale of the solutions. Near the surface the anelastic 
solution differs from the exact solution because of the dif- 
ferent boundary condition: Exact solutions verify SP = 0, 
a condition which transforms into £ = or u = for the 
approximate solution. 

2.3 Comments 

The foregoing developments show that the anelastic approx- 
imation applies under rather more general conditions than 
the subseismic approximation; indeed, the only requirement 



is the smallness of the Brunt- Vaisala frequency or, in other 
words, a large separation between the acoustic spectrum and 
the gravity spectrum. As this latter condition is often met 
in stars we can expect that the anelastic approximation per- 
forms better when applied to star models. 

Concerning the subseismic approximation, it is clear 
that it can be applied in the central region of a star but 
that surface layers should be avoided. Smeyers introduces 
the notion of boundary layer to describe the regions where 
his asymptotic solutions are valid. However, these bound- 
ary layers are somehow special since their thickness can be 
comparable to the radius of the star (in polytropes for in- 
stance) 1 . 

Broadly speaking, it turns out that the subseismic ap- 
proximation has a rather local character while the anelas- 
tic approximation has a global one. As eigenvalue problem 
are global problems in nature, the anelastic approximation 
should be better suited for these problems. 



3 EXAMPLES 

As in paper I we consider two polytropes: one of constant 
density and one of index n = 3. 



3.1 The homogeneous star model 

In this case analytic solutions exist either for the exact or 
the approximate equations (see paper I). 

In the asymptotic case of large wavenumbers (k — > oo), 
one finds that 



1 Classical boundary layers have a thickness very small compared 
to the size of the domain and which tends to zero as the small 
parameter is decreased. 
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Figure 2. Normalized eigenfunctions £ with £ = 2 and k = 10 (left) and fc = 20 (right) for the polytrope n = 3. As in Fig. 1, exact 
solutions are denoted by solid lines (with lo^ ~ 1.98 X 10~ 3 and oj| ~ 5.99 X 10 — 4 ) while the dashed and dotted lines correspond to 
their subseismic and anelastic approximations, respectively (with <^ 2 uba ~ 2.08 X 10~ 3 , w 2 nol ~ 2.02 X 10~ 3 for fc = 10 and c<j 2 ubs ~ 
6.15 x 10- 4 , ~ 6.04 x 10~ 4 for fc = 20). 



V^ 2 + > 



with A = 7 [k (f. + k + |) + £ + |] -2. Therefore 
2 2£(£ + 1) 1 



7 



2fc(2fc + 21 + 5) + 4£ + 6 - 4/7 ' 



Now using (22) of paper I we find that for the anelastic 
approximation 



2£{£ + 1) 



7 2k{2k + 2£ + 5) +4£ + 6' 
while the subseismic expression (21) can be rewritten as: 
i 2£(£ + 1) 1 



7 2fc(2fc + 3£ + 4 + 2/7)+6£+(2f + 4)/7' 

It is clear from these three expressions that, for high or- 
der modes, the anelastic approximation is very close to the 
exact expression. In fact, one finds that (lu 2 — wf ne i) /lo 2 ~ 
fc -2 while (lj 2 — cj 2 ubs )/o; 2 ~ k~ x , that is, the anelastic ex- 
pression converges quadratically while the subseismic one 
only linearly. 

This better behaviour of the anelastic approximation 
is confirmed by the shape of the eigenfunctions as shown in 
figure 1. There, we clearly see that the subseismic solution is 
good only in the central regions (r < 0.2) while the anelastic 
approximation remains close to the exact solution almost to 
the surface. 



3.2 The polytrope n — 3 

For a polytrope n — 3 a similar behaviour exists although 
the difference between the two approximations is less pro- 
nounced. In I we observed that the eigenfrequencies con- 
verged at different rates, the anelastic approximation one 
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Figure 3. A schematic picture of the modes of a star viewed 
from the anelastic viewpoint (above) and subseismic viewpoint 
(below), wjv and ui c are respectively the frequency of the lowest 
order gravity and acoustic modes. 



converging faster. Here we plot two eigenmodes of high order 
(fc = 10 and fc = 20) computed with the two approximations 
and with the complete equations (see figure 2). As expected, 
while both approximations describe the central regions very 
accurately, the anelastic one remains closer to the 'exact' 
solutions on a larger volume. For the fc = 20-mode, it de- 
parts noticeably from the exact solution close to the surface 
(r ~ 0.83). 



4 CONCLUSION 

In this paper we tried to clarify the differences between the 
subseismic and anelastic approximations which both aim 
at describing the low frequency spectrum. The subseismic 
approximation appears when one concentrates on the low 
frequency high radial order modes in the central region of 
a star; no constraint is imposed to the Brunt- Vaisala fre- 
quency. 

On the other hand the anelastic approximation assumes 
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a weak stratification but imposes no constraint on the degree 
of the mode. 

Hence, while the anelastic approximation makes the 
Brunt- Vaisala frequency, and thus the frequency of all grav- 
ity modes, vanishingly small compared to acoustic frequen- 
cies, the subseismic approximation focuses on gravity modes 
whose radial order is very large and hence have small fre- 
quencies compared to acoustic ones. 

In other words, the anelastic approximation removes the 
elasticity of the fluid by rejecting acoustic frequencies to in- 
finity and therefore allows for a description of the full spec- 
trum of gravity modes while the subseismic approximation, 
keeping ui c and u>n in a finite ratio, concentrates on one part 
of the spectrum, namely that containing high radial order 
modes which are the least sensitive to the elasticity of the 
fluid. This situation is summarized in figure 3. 

Since in stars the situation is often that ujn <C ui c , the 
use of the anelastic approximation is recommended as it is 
likely closer to the solutions of the complete equations; on 
the other hand, the subseismic approximation may be useful 
when one needs an analytic expression of gravity modes in 
the central regions of a star. 

Finally, it is worth mentioning the work of Durran 
(1989) who discussed these two approximations in the con- 
text of atmospheric sciences. In this field, where the sub- 
seismic approximation is called the "pseudo-incompressible 
appoximation" and the anelastic approximation the "modi- 
fied anelastic approximation" , the subseismic approximation 
appears to be superior to the anelastic approximation as it 
conserves the energy, a property which is important for non- 
linear problems. This result shows that the best choice for 
filtering out acoustic modes is dependent on the problem 
at hands. Therefore, our results which favour the anelastic 
approximation when searching for low-frequency modes of 
stars, may be specific to eigenvalue problems. 
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